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SUMMARY 


To this point in time flight simulation used in flight training and 
research has employed ideal straight-line paths for the glide slope and 
localizer portions of the ILS. This work reported gives, instead, practical 
path shapes and structures for simulator use. The difficulty is that real 
ILS courses are not smooth paths but do possess noise which is presented to 
the pilot. The work reported here will improve the quality of flight simu- 
lation by supplying the pilot with a more realistic indication of ILS 
paths. 


Records of real-world ILS facilities have been obtained or made and 
these form the basis of the research product, viz, digitized real-world 
paths prepared for inclusion in the programs which the flight simulator 
uses for producing the indications on -the pilot's course deviation indica- 
tor (CDI). In addition, some synthesized paths derived from some FAA sup- 
ported parallel work [1] are provided to cause the pilot to experience path 
excursions at Category II tolerance limits. 

The computer programs and data presented in this report provide a 
flight training research capability not previously available. With the 
emphasis on the workload of a single-pilot IFR (SPIFR) situation, the abil- 
ity to simulate realistically the aircraft cockpit environment for investi- 
gating pilot performance is important. This work will allow more realistic 
flight simulations of ILS approaches. Typically in the past a flight 
simulator used perfectly smooth ILS paths as input to the graphic display 
device usually course deviation indicators. 



II. INTRODUCTION 


The idea of using practical, real-world ILS beam structures in flight 
simulators is believed to have considerable merit, in particular with 
respect to flight training. Investigation has revealed that while many 
pilots using a simulator believe practical, representative beam structures 
were used, some who have flown actual ILS Category III approaches are well 
aware that the simulator does not represent actual conditions. Recognizing 
this problem, NASA Langley Research Center established contract NAS1-17368 
with the Avionics Engineering Center of Ohio University to implement 
realistic ILS course structures for an aircraft simulator at Langley. The 
major program steps are: 

1. Determine where and to what extent other contractors may have 
prepared non— ideal structures for use in simulation. 

2. Review existing measured data from various ILS facilities and 
select ten sets of data for implementation in the simulator. 

3. Digitize the selected paths with the structure noise included. 

4. Investigate a means of implementing the roughness to provide 
the most efficient method and to produce the most realistic 
results in terms of the pilot seeing what actually exists 

in space at the particular facility. 

5. Prepare several generic paths to allow a rigorous approach 
to determine pilot difficulties. 

6. Subsequent to selecting a method of inputting the noisy path, 
prepare ten (10) actual paths in a format compatible to the 
NASA Cyber computer. 

7. Prepare a setup of a cockpit display at the contractor’s 
facility to allow inspection of the results of the implementation 
of the path with noise. 

Actions taken in the completion of these tasks are discussed in this report. 

Also, as a part of this work, inquiries were made of several large 
commercial manufacturers of aircraft simulators regarding ILS course struc- 
tures produced in simulators. Of the 5 companies contacted, none indicated 
that actual path structures were being used. The typical approach has 
been to assume a perfectly straight structure oriented properly with 
respect to the runway. 

In addition to the inquiries to industry, a literature search was 
conducted. Dialog (an electronic data base operated by Lockhead 
Information Services) was searched for relevant publications and papers. 
Nothing relevant was found. Further, the Avionics Engineering Center's 
library contains no reports of significant work in this area. 
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III. REVIEW OF EXISTING MEASURED DATA FROM ILS FACILITIES 


When precision measurements are made of localizer and glide slope 
structures, it is necessary to provide a reference against which the path 
is compared. Two different reference systems are in common use. The FAA 
uses an inertial platform which is updated by DME information and manual 
marks when the flight passes over the threshold. Ohio University employs a 
radio telemetering theodolite. This is an optical instrument which is used 
by skilled operator to track the airplane from a specially designated point 
on the ground. The angular position is taken in electrical analog form 
from the theodolite and telemetered to the aircraft. Here it is subtracted 
from the course indication produced by the ILS receiver giving a difference 
which is desirably independent of the aircraft position. 

By using the inertial reference or theodolite it is possible to 
remove for the most part the course information dependence on the aircraft 
position. Obviously, an aircraft flying above the course will get a fly- 
down indication. Unless the precise position of the aircraft is known, one 
cannot ascertain whether the aircraft was high or the path had a dip in it. 
The references eliminate this problem. 

It should be noted that the ideal glide slope on-course is a three- 
dimensional geometrical shape, which in the case of the common image system 
is a cone. If the shape of this cone is perturbed, typically the pertur- 
bation is not linear vertically and horizontally. Consequently, flight 
measurements made with the aircraft not on the on-course means that the 
actual path is predicted based on assumptions of linearity. For example, 
if the aircraft is off-course on the localizer for one measurement of the 
glide slope and for another the aircraft is on the localizer, then it is 
possible that different glide slope structures will be recorded. The 
reason is that the horizontal structure of the cone is not uniform. 

From these circumstances one can observe variations in path struc- 
tures that are even measured in close time proximity. This does not mean, 
of course, the structures are not repeatable, but rather that the flight 
tracks for measurements were not executed with sufficient precision. 

Precise positioning of the aircraft and use of an accurate reference system 
will insure repeatability in glide slope structure measurements. 

Data used in preparing the tracks for these simulations were scruti- 
nized for repeatability which confirms .use of a satisfactory reference and 
reasonable positioning of the aircraft while making measurements. 

Approximately AO FAA flight check recordings were obtained through a 
visit to the Flight Inspection Field Office based at Oklahoma City, 

Oklahoma on July 15, 1983. These localizer course structure recordings 
were made of many different ILS facilities in the U.S. The localizer 
course structure recordings were acquired through the use of the FAA's 
Automated Flight Inspection System or AFIS (see figure 1).[2] The AFIS 
uses position information supplied by an inertial navigation system in a 
comparison with the position indicated by the measured ILS RF fields. The 
result of this comparison is the difference between the actual and 



AFIS Instrument Landing System 



Figure 1. Localizer Course Structure Data Taken by AFIS, 


indicated position of the aircraft. This information represents the error in 
the localizer course structure. The recordings were carefully reviewed and 
the best ten were selected as a basis for the real localizer path error 
structures to be prepared for NASA. The paths were selected based on the 
significance of abnormalities displayed by each path. Those paths with the 
greatest deviations being selected. 

Glide slope courses were selected from ILS flight test measurements 
made by the Ohio University Avionics Engineering Center staff between 1980 
and 1982 (See figure 2). Course structure measurements are accomplished 
by subtracting the known position of the aircraft from the observed (CDI) 
indication. The aircraft's position is tracked accurately using a radio 
telemetering theodolite. This information is then fed, along with the 
aircraft CDI indication, to a differential amplifier. The output of this 
amplifier is the difference between the known position and indicated posi- 
tion of the aircraft. This is the glide slope structure error information 
to be utilized by the simulation. Approximately 100 recordings were 
reviewed. Fifteen of the most interesting glide slope structures were 
selected and digitized for further inspection. Ten courses of these fif- 
teen were selected for use. 

A listing of the courses selected and the installations from which 
the data was taken is given in appendix A. Plots of the localizer and 
glide slope courses are also given in appendix B. Table B-I, located in 
appendix B, specifies which glide slope and localizer courses were paired 
together to form the ILS courses. 
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Figure 2. Glide Slope Course Structure Data Taken by the Staff of Avionics Engineering Center 























































IV. DIGITIZING THE SELECTED PATHS WITH STRUCTURE NOISE 

The courses selected were digitized using an analog strip-chart to 
digital data translator developed by the Avionics Engineering Center in 
1976. [3] The position of a manually-operated pen is sensed by the data 
translator and is coded into BCD format for recording on a 9-track digital 
magnetic tape recorder. Sample rates of the data translator are adjustable 
to run the manual track at a comfortable speed, allowing good resolution 
s and accuracy in the converted data. Sensitivity of the device is also 

adjustable to allow for the most efficient use of the available quanti- 
zation. The best settings for both the sample rate and the sensitivity 
were determined empirically and the same settings were used while digi- 
tizing all courses. 

The tape with the courses encoded into a BCD format was delivered to 
Ohio University Computer Services. The data were transferred to an IBM 370 
mainframe and stored for further processing. It was necessary to develop 
software to read the tapes and to unpack and translate the BCD coded data. 

A listing of the software used can be found in appendix C. The data for 
each localizer course were scaled to represent degrees deviation from the 
runway centerline. The software used in the rescaling process took into 
account the scaling of the original measurements, the sensitivity setting 
of the data translator and the localizer course width of each particular 
site. All ten courses were scaled such that the error will have the same 
relative significance on the four-degree course width (assumed by the NASA 
simulator) as it did with respect to the course width of the original ILS 
site. Glide slope course structures were scaled similarly. The NASA simu- 
lator assumes a 0.7 degree deviation from the glide slope angle for full 
scale deflection qf the CDI needle. The glide slope data were scaled to 
account for this as well as the scaling of the original measurements and 
the data translator used to digitize the data. 

An interpolation routine was developed which would readjust the 
sampling of each course structure such that samples occurred at 25-foot 
intervals. This normalization was applied to all ten courses so that all 
ten would be uniform in this respect. Specifying a uniform course format 
simplifies the implementation of the course stuctures in the existing NASA 
software, as the routine will be the same regardless of which particular 
course structure is used. The original measurements were made at various 
chart recorder speeds. Event marks corresponding to nautical mile marks 
were recorded during the data translation and this information was -used to 
readjust the courses into the standardized form. 

Each course was truncated to 1000 points (hence 25000 feet of local- 
izer and glide slope structure or just over 4 nautical miles). The last 
500 feet of course data (farthest from the runway) was altered so that the 
error data for both localizer and glide slope rises smoothly from zero to 
its actual value at 24,500 feet from the runway threshold. This was done 
to avoid a discontinuity in the CDI indication when the simulated aircraft 
first enters the area within 25,000 feet of the runway threshold for which 
■ the course structure data exists. The error beyond 25,000 feet is assumed 

to be insignificant for this study. 
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V. 


IMPLEMENTING THE ROUGHNESS 


The noise structure of ILS courses in space, whether localizer or 
glide slope are repeatable when measured precisely. The structure contains 
noise which is a result of the multipath effects of objects surrounding the 
airport such as buildings, terrain, power lines and towers. An aircraft 
moving along a given path in space will, provided no changes in terrain or 
building or structure changes have occurred, receive the same structure 
information. In fact, this characteristic quality of each localizer and 
glideslope is used to ascertain that if on successive fight evaluations 
the characteristic structure of the path is the same then the flight 
measurement is valid. The structures presented in this report are true 
path structures as would be seen in actual flight. They are not due to 
noise in the receiver telemetry or recording systems. These noise levels 
are below the observable thresholds. 

Several schemes for implementing the course roughness were discussed. 
The simplest method was determined to be a direct table lookup. The 
distance of the aircraft from the runway threshold was used in a formula to 
calculate the table index. This method is economical in that it is memory- 
conservative as well as easily implemented. For this table lookup approach 
it is necessary that all the courses assume a standard form. Data tables 
of 2,000 elements, each representing the localizer and glide slope struc- 
tures sampled at 25-foot intervals, were created. The tables are arranged 
as 1,000 pairs of elements, where each pair corresponds to the localizer 
and glide slope error in degrees for a particular distance from the runway 
threshold. Distance is determined by the position of the pair in the 
table. Figures 3 and A illustrate the relationship between the values in 
the table and the distance from the runway threshold. The model assumes 
bilateral uniformity. At any given point it is assumed that the error in 
the CD I indication is a function of the aircraft's distance from the runway 
threshold and not a function of the aircraft's lateral position with 
respect to the runway centerline. 

The table for a specific site may be loaded prior to the beginning of 
the simulation. This information is then referenced throughout the flight. 
For distances beyond 25000 feet beyond the runway threshold, the last ele- 
ment in the table will be selected (which is zero error, as described 
above). As the distance to the threshold falls below 25,000 feet the table 
lookup routine will index farther down the table until, as the aircraft 
crosses the threshold, the first pair of elements in the table will be 
selected. The values resulting from this table lookup process are simply 
added to the localizer and glide slope indication which is calculated by 
existing NASA simulator software. The calculated ILS CDI indication is, of 
course, a function of the aircraft lateral displacement with respect to the 
runway centerline and displacement with respect to nominal glide slope 
angle. However, the ILS course structure error information is only a func- 
tion of the distance to the runway threshold. 

For each of the various ILS simulations, approximately four miles of 
localizer and glide slope data are contained in a table of 1 ,000 pairs of 
elements. Each pair of elements corresponds to the error in the Indicated 
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position in degrees at intervals of 25 feet. The table is constructed such 
that localizer data appears first. Hence, for n in the range 1 to 1,000, 
the data contained in the (2*n)-l table position corresponds to the error 
in the indicated position (in degrees from centerline) at a distance of 
(n-l)*25.0 feet from the runway threshold for the localizer. Similarly, 
the 2*n element in the array corresponds to the error in the indicated 
position in degrees from the glide slope reference angle at a distance of 
(n-l)*25 ft from the runway threshold. The model for the course irregu- 
larities assumes lateral uniformity. The localizer values in the table 
correspond to the indicated position of an observer standing on the runway 
centerline at a distance x=(n-l)*25.0 from the runway threshold. 

The glide slope values in the table correspond to the indicated 
postion of an observer on the nominal glide path at some distance 
x=(n-l)*25 *0 feet from the runway threshold. Lateral uniformity implies 
that the relative error in indicated position at a given distance x will be 
the same for any position along the y axis at that distance (within the 
relevant range of y values). Since the maximum angular deviations are 
small this should be a valid approximation. The assumption of lateral uni- 
formity allows simple use of the x distance from the threshold to calculate 
the table index. (The n referred to above is the table element pair 
index.) 

Figure 5 is a section of a listing of the existing NASA software. 

Some modification to the subroutine RADNAV will be required. This section 
of the subroutine is shown below. In the existing software the variable 
XV0R1 corresponds to the distance along the x-axis from the runway 
threshold to the current position of the aircraft. This distance along 
with YV0R1 (the distance along the y-axis from the runway centerline) is 
used to calculate the CDI indication in degrees for the localizer. XV0R1 
is also used along with aircraft height above ground to calculate the glide 
slope indication. .The calculated glide slope deviation is referenced to 
the nominal glide path. 

To implement the non-ideal course structures it would be sufficient 
to alter the code shown in figure XX as follows: 


CURS CL= 1.0 
GSSCL-1.0 

C Specification of scale factors. This should be done during the program 
C initialization. 


C In the RADNAV section beginning at line 494: 


. 

L0C1=ATAN2(— YVOR1 ,ABS(XV0R1- TABLE ( 105)) )*RADDEG RADNAV494 
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C****** CALCULATE COMPASS MAGNETIC ERROR 
C 

CERR = (8.2103E-07 * X + 3.2994 5E-06 * 

I F ( LD I S I (34 ) ) CERR = 0. 

C 

IF (.NOT. NAVCOM2) GO TO 10 

COMPUTE V0R1 COORDINATES 

CALL STAT10NC ISTATN, NAVFRQ1 ,FREQ,NMAX) 
IF (ISTATN, GT. NMAX) GO TO 10 


RADNAV 442 
RADNAY 443 

Y + 7.965/57.295 RADNAY 444 

RADNAY 445 
RADNAV 446 
RADNAY 447 
RADNAV 448 
RADNAY 449 
RADNAV 450 
RADNAV 451 
RADNAV 452 


LOCIT = ICODE( ISTATN) .EQ. ILL 
Y0R1 = .NOT. LOCIT 

B0BS1 = ATAN2 (BV0R1S, BV0R1C) 

IF (.N0T.L0C1T)G0 TO 30 
BV0R1S = SIN(RUNWAY ( I ST ATN) 757.295 
BV0R1C = COS (RUNWAY ( I STATN5/57.295 
GO TO 31 
30 CONTINUE 


RADNAV 457 
RADNAV 458 
JFIX 19 
RADNAV 459 
RADNAV 460 
RADNAV 461 
RADNAV 462 
RADNAV 463 


SR = SI N(ORI ENT( I STATN )/57 .295) RADNAV 

CR = COS ( OR I ENT ( I STATN ) /57 .295 ) RADNAV 

BS = BV0R1S RADNAY 

BC =* BV0R1C RADNAV 

BV0R1S = BS*CR - BC*SR RADNAV 

BV0R1C = BC*CR + BS*SR RADNAV 

31 CONTINUE RADNAV 

DELX1 = X - XSTN( ISTATN) RADNAV 

I FCDELX1 .EQ 0.) DELX1 = 1.E-6 RADNAV 

DELY1 = Y - YSTN(ISTATN) RADNAV 

BITOSTN = ATAN2(-DELY1 ,-DELX1 ) RADNAV 

I F ( .N0T.L0CIT)B1T0STN=31T0STN+0RIENT( I STATN5/57.295 RADNAV 

RMAX = 1.E30 RADNAY 

IF (LOCIT) RMAX = 1.E30 - RADNAV 

GSRANGE = SQRT(DELX1 *DELX1+0ELY1 *DELY1 ) .LT. RMAX RADNAV 

XV0R1 = DELX1 *C0SV0R1 + DELY1*$I NV0R1 RADNAV 

YV0R1 = -DELX1 *$ INV0R1 + DELY1 *C0SV0R1 RADNAV 

TEMP = XVORI *XV0R1 + YV0R1*YV0R1 RADNAV 


DME1 = SQRT(TEMP + (H-ZSTN( ISTATN*(H-ZSTN( ISTATN) ) ) *0.000164468 RADNAV 


464 

465 

466 

467 

468 

469 

470 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 


DM E RATE = (DEMOLD - DME1 ) *3600. /H7 
DMERAT 1 = ABS (DMERATE) 

DMEOLD = DM El 
RANGE 1 = SQRT(TEMP) 

RCONE = 1 .428 1*CH-ZSTN( ISTATN) 

IF (LOCIT) 1,2 
1 IF (GSRANGE) 3 ,4 

3 LOCI = ATAN2 ( -YVOR1 , ABS CXV0R1 -TABLE ( 1 05 ) ) *RADDEG 
GS1 = XVORI .LT. 1000. 

GSDEV1 = 0.0 

IF(GSI) GSDEV1 * ATAN2 <H-ZSTN( I STATN) , -XVORI )*RADDEG -GLSI 
T01 = FR0M1 = .F. 

GO TO 5 


RADNAV 487 
RADNAV 488 
RADNAY 489 
RADNAV 490 
RADNAV 491 
RADNAY 492 
RADNAV 493 
RADNAV 494 
RADNAV 495 
RADNAV 496 
RADNAV 497 
RADNAY 498 
RADNAV 499 


Figure 5. Section of Listing Received from NASA Programmer. 

In this section of the program the relative distance to 
runway threshold and the CDI indication is calculated. 



C Calculation of localizer indication in degrees. 

GS1=XV0R1 .LT.1000 • 

GSDEV 1=0.0 RADNAV495 

IF(GSl) GSDEV = ATAN2(H-ZSTN( ISTATN) ,-XVORl )*RADDEG-GLSI RADNAV496 
C Calculation of glide slope. 

C The above lines are unchanged. 

LDEX=INT(ABS(XVOR/25 .0 ))*2-l RADNAV497 

C Calculation of the the lookup table index as a function of 
C the x distance from the runway threshold. 

IF (LDEX.GT. 1999) THEN LDEX=1999 RADNAV498 

IF (LDEX.LE.O) THEN LDEX=1 RADNAV499 

C If LDEX is greater than 1999 the aircraft is farther than 

C 25000 feet from the runway threshold. If this is the case, the 

C last element in the table is selected. In all the tables this 

C last element is zero and therefore no error is added for 

C distances beyond 25000 feet from the threshold. 

L0C1=L0C1+ILSERR( LDEX)*CURSCL RADNAV500 

C The localizer indication is simply the sum of the calculated 

C localizer indication and the error data. • 

C similarly for the glide slope. 

GSDEV1=GSDEV1+ILSERR(LDEX+1)*GSSCL RADNAV50 1 


Where: ILSERR is an array of 2000 elements which is loaded 

with the localizer and glide slope error structure 
prior to the beginning of the simulation. 

CURSCL and GSSCL are optional variables which will 
allow scaling of the sensitivity of the system to 
the course width. If CURSCL=1.0 the error added 
to the calculated Indication is exactly the error 
as measured at the various sites. If for example 
CURSCL=1.01 then the error added to the calculated 
localizer reading will be 1% greater than measured 
and so on. This should allow some study of 
tolerable course irregularities. These variables 
should be initialized to a value in some convenient 
manner earlier in the program execution. 

The above approach assumes that the course structures can be loaded 
into a one-dimensional array of 2,000 elements (from tape or other storage 
media) prior to the beginning of the simulation run. All the data 
collected and digitized have been rescaled so that simulated errors have 
the same relative significance on the 4-degree course width assumed by the 
NASA simulator for the localizer and the -K).7 degree angle (referenced to 
the nominal glide path) assumed for the glide slope. Figure 6 is a section 
of a listing of the existing NASA software. This section of the code out- 
puts the calculated CDI indication to the appropriate D/A converter to 
drive the instruments in the cockpit. In this section the data are scaled 
such that full scale deflection occurs for +2 degrees deviation from cen- 
terline for the localizer and -+0.7 degrees deviation from the nominal glide 
path for the glide slope. The variable CURSCL would allow the sensitivity 
of the real localizer signal to be altered if desired. In a similar manner 
GSSCL would allow alteration of the glide slope sensitivity. 
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GLIDE SLOPE FOR HSI ON 

DACOUT 419 

C 


DACOUT 420 


DACC30) = -GSDEV1 * SDAC19 

DACOUT 421 

c 


DACOUT 422 

c** 

HSI VOR/LOC INDICATOR 

DACOUT 423 

c 


DACOUT 424 


DACC28) = LOCI *SDAC16 

DACOUT 425 


1 F ( (LOCIT) ) DACC28) = ,4*L0C1 

DACOUT 426 

c 


DACOUT 427 

c 


DACOUT 428 

c* 


DACOUT 429 

Q********************** 

DACOUT 430 

c ********************* signal for sound system 

DACOUT 431 


Figure 6. Section of Listing Received From NASA Programmer. 
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VI 


GENERIC PATHS 


In addition to the ten real courses digitized for use in the simula- 
tor, nine generic rough paths were developed. Three sets of generic local- 
izer and glide slope paths were designed by Dr. R. H. McFarland, director 
of Avionics Engineering Center. The nine generic courses were created by 
taking all possible combinations of the three different idealized localizer 
courses with the three different glide slope courses. The generic local- 
izer and glide slope courses were designed to provide interesting noise 
disturbances for the pilot/aircraft control loop. These courses will be 
useful in a study of specific types of pilot difficulties when flying 
instrument approaches in non ideal situations. The generic courses were 
combined in the following manner: (See the plots accompanying this 

document. ) 
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Plots of the courses GENLOCn and GENGSn are given in appendix B. 
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VII, COCKPIT DISPLAY 


To facilitate assessment of the course structures chosen and 
designed, a cockpit display was fabricated in the laboratory. A flight 
simulation program was used in conjuction with a routine which simulates 
the table lookup technique discussed earlier in this paper to access the 
course structure data and simulate in-flight conditions. The flight simula- 
tor was directed to fly a perfect path (straight down the runway center- 
line, on the nominal glide slope angle). During the flight, the instrument 
indication was driven using the table lookup method. This allowed the 
error data (and only the error data) to be viewed. The results of these 
simulations were then downloaded to a Heath H-89 microcomputer. A simple 
routine was written in FORTH to take the data and output it at an appro- 
priate rate to a serial digital-to-analog converter. This serial device 
has two channels of D/A conversion available. [4] These D/A channels were 
used to drive a modified CDI instrument. The data were scaled and con- 
verted to the appropriate form during the initial simulation run to give 
the displayed data the same significance it will have on the NASA simula- 
tor. All the courses prepared for this study were reviewed on this 
display. 
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VIII. DIGITIZATION AND COURSE CONSTRUCTION SOFTWARE 


A variety of programs were written in a number of different languages 
to record, translate and display the data for the simulator^ ILS error 
structures. A description of the programs developed and their uses is 
given below. Listings of these programs may be found in appendix C. 

To facilitate the study of the effect of the structure noise on 
flight performance a FORTRAN program called USRPOSA was developed, which 
would simulate arbitrary flight paths. This program uses flight instruc- 
tion data to generate appropriate position information. This flight posi- 
tion information is then used by another routine which simulates the table 
lookup to be used in the NASA simulation. A FORTRAN program called FLY was 
developed to take the position data generated by USRPOSA and apply the 
table lookup method to generate the CDI indication information in much the 
same way as it will be done in the NASA simulator# This information was 
then reviewed on the laboratory cockpit display# Several IBM 370/CMS EXEC 
routines were developed to manage the file manipulations for each of the 
flight path programs in this study. 

The strip chart data translator used to digitize the measurements 
made at the various ILS sites produces a magnetic tape with all the samples 
written in BCD format, with each digit in a different byte# A simple CMS 
EXEC routine called TAMMOVE was written to allow these data files to be 
read by the IBM 370# Once the files had been read, it was necessary to 
translate the files from this packed BCD format to FORTRAN compatible 
integer format# This was done by using the program MAS5# MASS handles the 
conversion by using two variables which have been declared equivalent# One 
variable is of the logical type (1 byte) and the other is standard FORTRAN 
integer type. The program reads the digitized data and unpacks the records 
one digit at a time# After 4 digits have been unpacked, these four digits 
are used to calculate the integer number corresponding the BCD (4-digit 
representation). The most significant digit is used only as an event mark. 
This event information is preserved at this point as it will be needed 
later for spatial scaling of the data. A CMS EXEC was developed to define 
the files and determine the length of the input file (as this changes from 
one course to another) and pass this information on to the program MAS5# 

Once the data were obtained in integer representation, the next step 
was to scale the data. For this purpose two Pascal programs called SCLVERT 
and GSSCLV were developed. These programs are identical with the exception 
that one has the proper scaling constants for the localizer data (SCLVERT) 
and the other is designed to scale the glide slope data (GSSCLV). This 
difference in the scaling procedure arises from the fact that two dif- 
ferent measurement techniques were used in the original data acquisition. 
The localizer data taken by the FAA’s Automated Flight Inspection System is 
already expressed in degrees. The data for the glide slope paths taken by 
the Ohio University is in terms of microamperes# The end result of both 
programs is a 1 ,000 point table in a standard format with all values 
expressed in degrees. The original data were recorded at various chart 
speeds and hence, during the digitizing, the paths were translated with 
different sampling intervals. The scaling programs use an interpolation 
routine to adjust the spatial sampling to a uniform 25 feet per sample for 
all the error structures. Event marks corresponding to intervals of nauti- 
cal miles which were recorded during the data translation are used in the 



interpolation routine for the spatial scaling* Uniform spatial scaling was 
chosen in order to simplify the implementation of the noisy paths in the 
NASA simulator, since the method of implementation would be independent of 
the specific structure used for a particular flight simulation. The 
scaling of the magnitudes in the tables is also uniform in the sense that 
each point is scaled to have the same significance with respect to course 
width assumed by the NASA simulator as it had at the measured site. The 
program SCLVERT prompts the user for the course width of the original site 
in order to scale the localizer data accordingly. The glide slope data is 
already in terms of microamperes and hence this information is not needed. 

A CMS EXEC routine to make file definitions etc. , was developed for both 
scaling programs. 

Once the data was properly scaled and in a uniform format, a magnetic 
tape was created containing all ten course error structures. This tape was 
then delivered to NASA Langley with a short description of the recommended 
implementation method and details pertaining to the tape and data format. 
This information is repeated below. 



IX 


TAPE FORMAT 


The magnetic tape supplied to NASA as part of this contract is in the 
following format: 

9-Track 
1600 BP I 
Lrecl 80 
Block 800 
Recfm FB 
EBCDIC 


The ten courses were written on the tape one element per record 
using a FORTRAN format: 


FORMAT(F10 .8) 


The table is arranged in pairs of elements with the localizer data 
appearing first: 


Localizer error 
Glide slope error 
Localizer error 
Glide slope error 


etc. 


Each pair represents the noise at some distance from the threshold. 

The tape contains 10 real ILS course structures and 9 generic courses 
as described above. There are no tape marks between each course and no 
tape label at the end of all 19 courses, i.e.: 


Course 1 
Course 2 
Course 3 



Course 10 
Generic Course 1 
Generic Course 2 


Generic Course 9 

Each course is 2,000 elements long. Hence there are 38,000 records on the 
tape. 
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XI, 


APPENDICES 


1 • Appendix A * 

The data for the real localizer courses were taken from the following 
locations : 


LOCI 

- 

CDG 

L0C2 

- 

FWH 

L0C3 

- 

IIK 

L0C4 

- 

DRT 

L0C5 

- 

TIK 

L0C6 

- 

VQE 

L0C7 

- 

ALW 

L0C8 

- 

TCB 

L0C9 

- 

SPT 

LOC10 

- 

SPT 


Houston, Texas, ILS runway 32, Jan. 4, 1982. 

Carswell AFB, Fort Worth, Texas. Runway 35, July 13, 1982. 
Tinker AFB Oklahoma City OK. Runway 35, January 19, 1981. 

Del Rio, Texas. Runway 13, June 9, 1982. 

Tinker AFB Oklahoma City, Oklahoma. Runway 35, Jan. 19, 1981 
Randalf AFB, San Antonio, Texas. ILS runway 32L 
Walla Walla, Washington. ILS runway 20 
Fort Worth, Texas. Runway 17, July 1, 1982. 

Albuquerque, New Mexico. ILS runway 8, July 26, 1982. 
Albuquerque, New Mexico. ILS runway 8, July 26, 1982. 


The real glide slope course data was obtained from measurements made at the 
following times and places: 


GSl 

- 

SHV 

GS2 

- 

SHV 

GS3 

- 

SHV 

GS4 

- 

SHV 

GS5 

- 

SHV 

GS6 

- 

IPT 

GS7 

- 

IPT 

GS8 

- 

SHV 

GS9 

- 

IPT 

GS10 

- 

SHV 


Shreveport, Louisiana. Runway 13, May 27, 1981. 

Shreveport, Louisiana. Runway 13, May 27, 1981. 

Shreveport, Louisiana. Runway 13, May 28, 1981. 

Shreveport, Louisiana. Runway 13, May 28, 1981. 

Shreveport, Louisiana. Runway 13, May 28, 1981. 
Williamsport, Pennsylvania. Runway 27, Oct. 30, 1980. 
Williamsport, Pennsylvania. Runway 27, Oct. 21, 1980. 
Shreveport, Louisiana. Runway 13, May 28, 1981. 
Williamsport, Pennsylvania. Runway 27, Oct. 30, 1980. 
Shreveport, Louisiana. Runway 13, may 28, 1981. 



2. Appendix B 
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Table B-1 . Glide Slope Localizer Courses Paired Together to Form 
ILS Courses. 


Coursel ........... LOCI 

Course 2 ........... L0C2 

Course3 L0C3 

Course4 L0C4 

Course5 L0C5 

Course6 L0C6 

Course7 L0C7 

Course8 L0C8 

Course9 L0C9 

CourselO LOCIO 


23 
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0.00 
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Figure B-l* Generic Glide Slope No* 1 (GENGS1)* 
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Figure B-4. Generic Localizer No. 1 (GENLOC1). 
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Figure B-17. Localizer Course 1 (LOCI) 
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Figure B— 18 . Localizer Course 2 (LOC2) 
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Figure B— 19. Localizer Course 3 (LOC3) 
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Figure B— 21 . Localizer Course 5 (LOC5) 
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Figure B-22. Localizer Course 6 (LOC6) 
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Localizer Course 8 (LOC8) 
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Figure B-25. Localizer Course 9 (LOC9) 
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Figure B-26. Localizer Course 10 (LOCIO) 


3. Appendix C. Program Listings 


This exec allows the IBM370 to read tapes created by the strip chart to digital 
data translator. 


FI INMOVE TAPI (RECFM U LRECL 20000 BLOCK 20000 
FI OUTMOVE DISK A A C (RECFM U LRECL 20000 BLOCK 20000 
MOVEFILE 
XEDIT A A C 
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4C0NTR0L OFF 

4IF .41 EQ .7 4GOTO -TELL 

LISTFILE 41 42 43 (L NOHEADER STACK 

4READ VARS 4A 4B 4C 4D 4E 4F 4G 4H 41 4J 4K 

FI 3 DISK 41 42 43 (RECFM V LRECL 20000 BLOCK 20000 

FI 7 DISK C41 42 43 

FI 9 TERM 

4 TYPE FILE LENGTH 4E 
4 STACK 4E 
EXEC LODE MAS 5 

4TYPE CONVERSION COMPLETE.. .FILE C41 42 43 CREATED 

4EXIT 

-TELL 

4BEGTYPE 

* ...THIS EXEC TAKES A RECORD AS PRODUCED BY THE ANALOG TO DIGITAL 

* DATA TRANSLATOR AND CONVERTS IT TO A MORE USEFUL FORM. 

* THE EXEC DETERMINES THE LENGTH OF THE FILE, SETS UP THE PROPER 

* FILE DEFINITIONS AND THEN CALLS THE PROGRAM MAS5. MAS5 CONVERTS 

* THE FILE FROM BCD TO INTEGER REPRESENTATION. THE EXEC PASSES 

* THE LENGTH OF THE VARIABLE LENGTH RECORD TO THE PROGRAM . 

* USE OF THIS EXEC CREATES A NEW FILE AND DOES NOT DESTROY THE 

* INPUT RECORD. EXEC ASSUMES 1 RECORD IN FILE. FILE CREATED HAS 

* SAME FN FT FM' EXCEPT A »C> IS ADDED TO THE BEGINNING OF THE FILE 

* NAME. 

4 END 
4EXIT 



C THIS PROGRAM CONVERTS FILES WRITTEN IN BCD FORMAT BY THE ANALOG TO 
C DIGITAL DATA TRANSLATOR TO FORTRAN-PASCAL COMPATIBLE INTEGER FORMAT. 
C THE PROGRAM REQUIRES THE LENGTHT OF THE RECORD TO BE PROCESSED TO 
C BE CONTAINED IN 14 FORMAT AT THE BEGINNING OF THE FIRST RECORD. 

C USE OF THE CONVERT EXEC IS RECOMMENDED. 

INTEGERS L (4000), EVENT, N3,N2,N1 
LOGICAL*1 K (2 0000) 

LOGICAL*! LBUFF(4 ) 

EQUIVALENCE (LBUFFC1 ) ,INTGER) 

C THE PACKED BCD FORMAT IS UNPACKED INTO INTEGER NUMBERS 1 BYTE 
C AT A TIME THROUGH THIS EQUVLENT LOGICAL/I NTEGER VARIABLE. 

READ (9 , 10)NN 

C LENGTH OF VARIBLE LENGTH RECORD IS READ 
10 FORMAT (15) 

READ (3 ,300 ,END=999) (K ( I ) , I =1 , NN) 

C ENTIRE RECORD IS READ INTO L0GICAL*1 ARRAY. 

300 FORMAT (1 50(1 0(1 1 A1 ) ) ) 

999 CONTINUE 
N=l/4 

C NUMBER OF INTEGER NUMBERS EXPECTED TO RESULT FROM BCD INPUT OF LENGTH 
C NN IS CALCULATED. 

WRITE (7 ,30)N 

C NUMBER OF VALUES IN THE OUTPUT FILE IS THE FIRST RECORD IN THE OUTPUT 
C FILE. 

DO 20 1=1 ,N 

C DO THE CONVERSTION FOR AS MANY VALUES AS EXPECTED. 

I NTGER=0 

C CLEAR OUT EQUIVALENT VARIABLE. 

I I = 1*4-3 

C INDEX BY FOUR AND SKIP THE FIRST BYTE. 

LBUFF(4)=K( 1 1+1 ) 

C STEP THROUGH BYTE BY BYTE FOR FOUR BYTES 
EVENT=I NTGER 
LBUFF(4)=K(I 1+2) 

N3=l NTGER 
LBUFF(4)=K(I 1+3) 

N2=INTGER 
LBUFF(4)=K( 1 1+4) 

N1=INTGER 

L( I )=EVENT*1000+N3*100+N2*10+N1 

C CALCULATE THE INTEGER EQUIVALENT OF THE LAST FOUR BYTES OF BCD 
C CODED DATA 
20 CONTINUE 

WRITE(7,30) (L( I ) ,1=1 ,N) 

C PLACE IN OUTPUT FILE 
30 FORMATO ',1016,/) 

1000 F0RMAT(1X,1 10) 

END 
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<« 

I***************************************************"********************** 

* * 

* THIS PROGRAM IS DESIGNED TO PROCESS COURSE STUCTURES FOR THE NASA * 

* SIMULATION PROJECT. THE PROGRAM IS DESIGNED TO READ FILES CREATED * 

* BY THE 'CONVERT* PROGRAM. THIS PROGRAM THEN TAKES THE VALUES AND * 

* SCALES THEM TO DEGREES (SEE DOCUMENTATION ON DIGITIZING PROCESS) * 

* THEN THIS PROGRAM USES THE EVENT MARKES IN THE CONVERTED FILE TO * 

* DETERMINE THE 01 STANCE OF EACH POINT FROM THE RUNWAY THRESHOLD. * 

* (IN FEET). NEXT THE PROGRAM CREATES A FILE OF 1000 POINTS WHICH * 

* THROUGH AN INTERPOLATION PROCESS HAVE BEEN ADJUSTED TO BE 25 FEET * 

* APART. THIS PROVIDES A COURSE STUCTURE WHICH IS JUST OVER 4 NM LONG. * 

* FINALLY THE VAULES IN THIS TABLE ARE AGAIN RESCALED BY: FIRST * 

* NORMALIZING BY THE WIDTH OF THE LOCALIZER PATH AT THE ACTUAL * 

* INSTALLATION FROM WHICH THE DIGITIZED DATA WAS TAKEN, AND THEN * 

* MILTIPLING BY THE LOCALIZER PATH WIDTH WHICH IS ASSUMED BY THE * 

* NASA SIMULATOR (4 DEGREES). THIS RESCALING ALLOWS THE DATA POINTS IN* 

* THE TABLE TO HAVE THE SAME (RELATIVE) SIGNIFICANCE AS THEY WOULD AT * 

* THE INSTALLATION. THE USER IS PROMPTED TO ENTER THE LOCALIZER PATH * 

* WIDTH FOR THIS SCALING OPERATION. * 

* * 
*#****#******#***m******#**#*»*m**mt*tm***mt**mmt#*»*m***) 

PROGRAM SCALVERTC INPUT f OUTPUT); 


CONST 

(* DX IS THE INCREMENTAL DISTANCE FOR THE INTERPOLATION *) 
DX=25.0; 


VAR 

NUM: I NTEGER; (* NUMBER OF SAMPLES IN IfPUT *) 

ECNT: INTEGER; (* NUMBER OF EVENT MARKS *) 

POSITION:ARRAY(.U. 5000,1. .2.) OF REAL; 

(* POSITION AND ANGLE STORAGE FOR FIRST PASS *) 
RAWSTUF :ARRAY (.1 ..5000.) OF REAL; 

(* RAW, UNSCALED INPUT DATA FROM CONVERT *) 
1,11,111:1 NTEGER; (* GLOBAL I NDEXERS *) 

EVENT : ARRAY (.1.. 30.) OF INTEGER; 

(* STORAGE FOR EVENT POSITIONS *) 

COR: ARRAY ( .1 . .1000.) • OF REAL; 

(* FINISHED TABLE STORAGE SPACE *) 
DPERSAM:REAL; (* DISTANCE PER SAMPLE ON INPUT *) 

( *-.*** ******************************* *+****■**+******* *************** 

* RAWIN. ...THIS PROCEDURE READS IN THE DATA FILE CREATED BY CONVERT * 

* AND STORES IT IN RAWSTUFF ARRAY. * 

**• *******************•»*****»********■»*******************************) 
PROCEDURE RAWIN; 

VAR 

1,11,111:1 NTEGER; 

EFLOP: BOOLEAN; 


BEGIN 

EFLOP :=FALSE; 

READ (NUM); 

ECNT:=0 ; 

FOR I :=1 TO NUM DO 
BEGIN 

IF NOT EOF (INPUT) THEN READ ( RAWSTUF (.1 .) ); 

(* READ RAW DATA I WTO RAWSTUFF ARRAY *) 

IF (RAWSTUF (•!•)>! 000) AND ( (RAWSTUF( .1 . )-1 000) >1000) THEN RAWSTUF(.U) 
RAWSTUF ( *1-1 •); 

(♦ IF INCOMMING DATA I S OUT OF RANGE AND IT IS NOT AN EVENT MARK, 
THEN IGNORE AND SET CURRENT SAMPLE EQUAL TO LAST SAMPLE *) 

IF RAWSTUF ( • l • )>1 000 THEN 

(* IF IT MAKES IT HERE IT MUST BE AN EVENT *) 


(* INDEXES *) 

(* EVENT FLIP FLOP, USED TO INSURE 
EACH EVENT IS COUNTED ONLY ONCE *) 


(* NO ACTIVE EVENT *) 

(* READ NUMBER OF SAMPLES *) 
(* BEGIN COUNTING EVENTS *) 
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END; 


BEGIN 

<* PROCESS THE EVENT MARK *) 

RAW STUF (•!.): =RAW$TUF ( . I . )- 1 000 ; 
IF EFLOP=FALSE THEN 
BEGIN 

ECNT :=ECNT+1 ; 

EVENT < .ECNT. ) : = l ; 
EFLOP:=TRUE; 

END; 

END 


ELSE EFLOP:=FALSE; 

(* DONE IN T HE EVENT OF NO EVENT.... 
RAWSTUF ( . I • ) :=RAh. ( . I . >-500; 

(* SUBTRACT OUi DC OFFSET *) 

END; 


*) 


{•#*■»******»■»*****■**■*****■**«***********»»■»*■» 

* ANGLESCALE THIS PROCEDURE SCALES THE VALUES IN POSITION TO * 

* THE CORRESPONDING ANGULAR VALUE AS DETERMINED BY THE * 

* DIGITIZING METHOD. * 

4t4tMt « At « t « t **«Mtt****ft***ft**t***m***m********«********«***) 

PROCEDURE ANGLESCALE; 

YAR 1,1 1,1 1 1 : INTEGER; 

BEGIN 

FOR 1 5*1 TO NUM DO 
BEGIN 

POSITION( .1 *2.) :=RAWSTUF(. I .) *(-0.005); 

(* EACH BIT OF QUANTIZATION CORRESPONDS TO 0.005 DEG. *) 

END; 

END; 


* DISTANCE.. ..THIS PROCEDURE USES THE EVENT MARKES (COUNTED AND * 

* LOCATED IN RAW IN PRXEDURE) TO DETERMINE THE DISTANCE VALUES * 

* FOR POSITION ARRAY * 

*t1MHHHM(#<MhHHW..#*##t.**..*t*#.*1Hl************************************* ) 

PRXEDURE DISTANCE; 


VAR 

1,1 1,1 11:1 NTEGER; 

AVGSMP:REAL; (* AVERAGE NUM OF SAMPLES PER NM *) 

SUM .MAXPOS: REAL; (* SUM USED IN AVG PRXESSING 

MAXPOS.. POSITION OF FARTHEST SAMPLE *) 


BEGIN 

SUM:=0; 

FOR I :=2 TO ECNT DO 

SUM : = EVENT ( . I . )-EVENT( . I - 1 . )+SUM; 

AVGSMP:=SUM/(ECNT-1); C* AVG SAMPLES PER NM FOUND *) 

DPERSAM:=6076.0/AVGSNP; (* FEET PER NM FOUND *) 

MAXPOS : “EVENT ( .ECNT. )*DPER$AM; 

FOR I :«1 TO NUM DO 

POSITIONS 1 , 1 . ) :*MAXPOS-< 1-1 )*DPERSAM; 

(* CALX LATE POSITIONS WORKING BACK FROM THE MAXIMUM POS.*) 


END; 

(m#tmt**#**#t*m#***#mti»»********«******w*#t*t*#***#m# 

• STRUCTUREOUT.... THIS PROCEDURE WRITES OUT A FILE OF THE * 

* FIRST PASS RESULTS...^ 

PROCEDURE STRUCTUREOUT; 
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VAR 

1,1 1 : INTEGER; 


BEGIN 

WRITELNP * ,NUM:7); 

FOR I :=1 TO NUM DO 

WRITELNP '.POSITION!. I ,1.):10:2,' ' ,POSITION(.l ,2.) :10:3); 


( *******1MM>****tt****»****t**t*********t***t******t**************** 

* UN I COR THIS PROCEDURE CREATES A ’UNIFORM' TABLE BY WAY OF * 

AN INTERPOLATION ROUTINE. ALL TABLES HAVE 1000 ENTRIES * 

* WHICH ARE 25 FEET (DX) APART. * 

*****************■*•*»***********************************************) 


PROCEDURE UN I COR; 


VAR 

SL,XX,DXX,CURPOS:REAL; 

(* SL - SLOPE FOR STRAIGHT LINE INTERPOLATION APPROXIMATION *) 
I* XX - DISTANCE BETWEEN POINT TO BE INTERPOLATED AND THE 

LOWEST KNOWN BRACKETING VALUES *) 

<* DXX - I NCRENTAL Y VALUE FROM I NTERPOLATION *) 

(* CURPOS - I NSTANTAENOUS POSITION TO BE INTERPOLATED FOR *) 
BACK.NOPOI NT :BOOLEAN; 

(* BACK - FLAG FOR BACK COURSE *) 

(* NOPOINT - FLAG WHICH INDICATES THAT THE CURPOINT IS NOT 
BETWEEN THE BRACKETING VALUES *) 

1,11,111:1 NTEGER; 


BEGIN 

1 1 :=1 ; (* START COUNTING *) 

CURPOS: =25000.0; (* START FROM FARTHEST POSITION *) 

FOR I :=1 TO 1000 DO 
BEGIN 

N0P01 NT:=TRUE; 

WHILE NOPOINT DO 
BEGIN 

IF CURPOS > POSITION!. I 1,1.) 

THEN NOPOINT:=FALSE (* CURPOS BETWEEN B POINTS *) 

ELSE lls-ll+1; (* ELSE MOVE TO NEXT SET OF POINTS *) 
END; 

(* WHEN THIS POINT IS REACHED A PROPER POINT FOR INTERPOLATION 
HAS BEEN FOUND, SO BEGIN INTERPOLATING *) 

SL: = (POSITION(.l 1+1 ,2.)-P0SITI0N(. I 1 ,2.))/DPERSAM; 

XX:=CURPOS - POSITION!. I 1,1.); 

DXX:=XX*SL; 

I 1 1 :=1 001-1 ; 

COR! .1 11.) :=POS IT10N! . I I 
CURPOS: =2 5000. 0-1 *25.0; 

END; 

END; 

( ************************************************************** 

* CORSCALE. ...THIS IS THE FINAL SCALING PASS. THE USER IS * 

* PROMPTED FOR THE LOCALIZER PATH WIDTH FROM THE * 

* ACTUAL SITE. THE TABLE IS NORMALIZED BY THIS * 

* ANGLE AND THEN RESCALED BY THE WIDTH ASSUMED BY THE * 

* NASA SIMULATION (4 DEGREES) * 

***♦****■»**•*****»*»********■•****'»****•***********************) 
PROCEDURE CORSCALE; 

VAR 

1,1 I : INTEGER; 

OWIDTH:REAL; (* ORIGINAL LOCALIZER WIDTH *) 


(* THESE STATMENTS REDIFINE THE FILES FOR USER INTERACTION *) 
CLOSE! INPUT); 

SYSTEM (’ FI SYS IN TERM’ , I ); 
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CLOSE(OUTPUT); 

SYSTEMCFI SYSPRINT TERM* , I ); 

WRITELNC ENTER LOCALIZER COURSE WIDTH FOR THIS INSTALLATION.. ».' ); 
RESET(lfPUT); 

READLN(OWIDTH); 

FOR I :=1 TO 1000 DO 

COR( . I . ) : =COR( . I . )/0W I DTH*4 .0; 


/ ***************** 

* UNICOROUT CREATES a file on DISK CONTAINING THE FINISHED 

TABLE. 

•**•••*« 

PRXEDURE UN I COROUT; 

VAR 

1,1 I -.INTEGER; 


BEGIN 

(* DEFINE FIDEFS FOR DISK FILE OUTPUT *) 
CLOSE (OUTPUT); 

SYSTEMCFI SYSPRINT DISK UC CHART C',tl); 
REWRITE(OUTPUT); 

FOR l:*1 TO 1000 DO 
WRITELN(COR(.l.)s10:8); 


* SMOOTHEND.....THIS PROCEDURE ’SMOOTHES* THE END OF THE TABLE * 

* SO THAT THE FIRST SMPLE WILL NOT CAUSE A BIG JUMP ON * 
CD I NEEDLE. 


PRXEDURE SMOOTHEND; 


VAR I : INTEGER; 
SLOPE :REAL; 


BEGIN 

SLOPE :**(-COR( .975.) )/625.0; 

FOR I :=976 TO 1000 DO 

(* draw a STRAIGHT LINE FROM THE 975TH SAMPLE TO ZERO *) 
BEGIN 

COR( .1 . ) :=XR( .975. )+( l-975)*25*SL0PE; 

END; 

END; 


(•MAIN PRXRAM BODY 


*) 


BEGIN 

RAW IN; 

ANGLESCALE; 

DISTANCE; 

STRUCTUREOUT; 

UNICOR; 

SMOOTHEND; 

CORSCALE; 

UNICOROUT; 

END. 
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<* 

********************** ********************************************* ****** 

* * 

* GLIDE SLOPE VERSION . * 

* THIS PROGRAM IS DESIGNED TO PROCESS COURSE STUCTURES FOR THE NASA * 

* SIMULATION PROJECT. THE PROGRAM IS DESIGNED TO READ FILES CREATED * 

* BY THE 'CONVERT' PROGRAM. THIS PROGRAM THEN TAKES THE VALUES AND * 

* SCALES THEM TO DEGREES (SEE DOCUMENTATION ON DIGITIZING PROCESS) * 

* THEN THIS PROGRAM USES THE EVENT MARKES IN THE CONVERTED FILE TO * 

* DETERMINE THE DISTANCE OF EACH POINT FROM THE RUNWAY THRESHOLD. * 

* (IN FEET). NEXT THE PROGRAM CREATES A FILE OF 1000 POINTS WHICH * 

* THROUGH AN INTERPOLATION PROCESS HAVE BEEN ADJUSTED TO BE 25 FEET * 

* APART. THIS PROVIDES A COURSE STUCTURE WHICH IS JUST OVER 4 NM LONG. * 

* FINALLY THE VAULES IN THIS TABLE ARE AGAIN RESCALED BY: FIRST * 

* NORMALIZING BY THE WIDTH OF THE GLIDE SLOPE PATH AT THE ACTUAL * 

* INSTALLATION FROM WHICH THE DIGITIZED DATA WAS TAKEN. * 

* THE USER IS NOT PROMPTED TO ENTER THE GLIDS SLOPE BEAM * 

* WIDTH FOR THIS SCALING OPERATION. * 

* * 

*************************************************************^**********^) 

PROGRAM SCALVERTC INPUT, OUTPUT); 


CONST 

(* DX IS THE INCREMENTAL DISTANCE FOR THE INTERPOLATION *) 
DX=25.0; 


NUM: INTEGER; (* NUMBER OF SAMPLES IN INPUT *) 

ECNT: INTEGER; (* NUMBER OF EVENT MARKS *) 

POSITION:ARRAY(.1.. 5000,1. .2.) OF REAL; 

(* POSITION AND ANGLE STORAGE FOR FIRST PASS *) 
RAWSTUF :ARRAY ( . 1 . .5000 • ) OF REAL; 

(* RAW, UNSCALED INPUT DATA FROM CONVERT *) 
1,11,111:1 NTEGER; (* GLOBAL I NDEXERS *) 

EVENT :ARRAY ( .1 ..30. ) OF INTEGER; 

(* STORAGE FOR EVENT POSITIONS *) 

COR: ARRAY (.1.. 1000.) OF REAL; 

(* FINISHED TABLE STORAGE SPACE *) 
DPERSAM:REAL; •. <* DISTANCE PER SAMPLE ON INPUT *) 


(******************************************************************* 

* RAWIN....THIS PROCEDURE READS IN THE DATA FILE CREATED BY CONVERT * 

* AND STORES IT IN RAWSTUFF ARRAY. * 

********************************************************************* 

PROCEDURE RAW IN; 


VAR 

1,11,111:1 NTEGER; 
EFLOP: BOOLEAN; 


(* INDEXES *) 

(* EVENT FLIP FLOP, USED TO INSURE 
EACH EVENT IS COUNTED ONLY ONCE *) 


BEGIN 


EFLOP :=FALSE; 
READ (NUM); 
ECNT:=1 ; 


(* NO ACTIVE EVENT *) 

(* READ NUMBER OF SAMPLES *) 
(* BEGIN COUNTING EVENTS *) 


FOR I :=1 TO NUM DO 
BEGIN 

IF NOT E0F( INPUT) THEN READ (RAWSTUF ( . I . ) ) ; 

(* READ RAW DATA INTO RAWSTUFF ARRAY *) 

IF (RAWSTUF( .1 ,)>1000) AND ( (RAWSTUF(.I .)-1000)>1000) THEN RAWSTUF(.I 


RAWSTUF (.1-1.); 

(* IF INCOMMING DATA IS OUT OF RANGE AND IT IS NOT AN EVENT MARK, 
THEN IGNORE AND SET CURRENT SAMPLE EQUAL TO LAST SAMPLE *) 
IF RAWSTUF (.1.) >1000 THEN 

(* IF IT MAKES IT HERE IT MUST BE AN EVENT *) 

BEGIN 

(* PROCESS THE EVENT MARK *) 



RAWSTUF ( . I . ) : “RAWSTUF ( . I . )- 1 000 ; 

IF EFL0P=FALSE THEN 
BEGIN 

EVENT(.ECNT. ):“l; 

ECNT:“ECNT+1 ; 

EFL0Ps=TRUE; 

END; 

END 

ELSE EFLOPs “FALSE; 

(* DONE IN THE EVENT OF NO EVENT.... *) 
RAWSTUF (.).): “RAWSTUF ( . I . )-500 ; 

(* SUBTRACT OUT DC OFFSET *) 

END; 

END; 


(********************«»***************•*************************** 

* ANGLESCALE THIS PROCEDURE SCALES THE VALUES IN POSITION TO * 

* THE CORRESPONDING ANGULAR VALUE AS DETERMINED BY THE * 

* DIGITIZING METHOD. * 

immuhmhhmhmmmhmmmmhhmmmmhmhmumm************************************) 

PROCEDURE ANGLESCALE; 

VAR 1,1 1,1 11 : INTEGER; 

BEGIN 

FOR I :«1 TO NUM DO 
BEGIN 

POSITION!. 1 ,2. ) :“RAWSTUF ( .1 . )*(-0.2); (* 2 MICO WP QUANTIZATION *) 

END; 

END; 


DISTANCE.... THIS PROCEDURE USES THE EVENT MARKES (COUNTED AND 

* LOCATED IN RAWIN PROCEDURE) TO DETERMINE THE DISTANCE VALUES * 

* FOR POSITION ARRAY * 


PROCEDURE DISTANCE; 


VAR 

1,1 1,1 ll:INTEGER; 
AVGSMP: REAL; 

SUM, MAXPOS: REAL; 


(• AVERAGE NUM OF SAMPLES PER W *) 

{* SUM USED IN AVG PROCESSING 

MAXPOS.. POSITION OF FARTHEST SAMPLE *) 


BEGIN 

SUM :=0; 

FOR I :=2 TO ECNT-1 DO 

SUM : “EVENT ( . I • )-EVENT ( . I -1 . )+SUM; 

AVGSMP: “SUM/CECNT-2); <* AVG SAtPLES PER NM FOUND *) 

DPERSAM:=6076.0/AVGSMP; (* FEET PER NM FOUND *) 

MAXPOS : =EV ENT ( .ECNT-1 . ) “OPERSAM; 

FOR I :“1 TO NUM DO 

POSITION! .1 , 1.) : “MAXPOS- < I -1 ) *DPERS AM ; 

(* CALCULATE POSITIONS WORKING BACK FROM THE MAXIMUM POS.*) 


END; 


( 

* 


* 

* 


STRUCTUREOUT.... THIS PROCEDURE WRITES OUT A FILE OF THE 
FIRST PASS RESULTS. 


PROCEDURE STRUCTUREOUT; 


VAR 

1,11 : INTEGER; 
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BEGIN 


WRITELNC 1 1 ,NUM:7); 

FOR !:»1 TO NUM DO 

WRITELNC 1 \POSITION(.t,U):tO:2 f , «, POSITIONS l ,2.) : 10:3); 

END; 

^ ****************************************************************** 

* UNICOR THIS PROCEDURE CREATES A 'UNIFORM' TABLE BY WAY OF * 

* AN INTERPOLATION ROUTINE. ALL TABLES HAVE 1000 ENTRIES * 

* WHICH ARE 25 FEET <DX) APART. * 

♦I.*****************************************************************) 

PROCEDURE UNICOR; 


VAR 


SL,XX,DXX,CURPOS:REAL; _ „ _ 

(* SL - SLOPE FOR STRAIGHT LINE INTERPOLATION APPROXIMATION *) 
(* XX - DISTANCE BETWEEN POINT TO BE INTERPOLATED AND THE 

LOWEST KNOWN BRACKETING VALUES *> 

(* DXX - INCRENTAL Y VALUE FROM INTERPOLATION *) 

(* CURPOS - I NSTANTAENOUS POSITION TO BE INTERPOLATED FOR *) 


BACK , NOPO I NT : BOOLEAN; 

(* BACK - FLAG FOR BACK COURSE *) 

{* NOPOINT - FLAG WHICH INDICATES THAT THE CURPOINT IS NOT 
BETWEEN THE BRACKETING VALUES *) 


1,11,111:1 NTEGER; 


BEGIN 

1 1 :=1 ; (* START COUNTING *) 

CURPOS :=25000.0; {* START ROM FARTHEST POSITION *) 

FOR I :=1 TO 1000 DO 
BEGIN 

NOPOINT:=TRUE; 

WHILE NOPOINT DO 
BEGIN 

IF CURPOS > POSITIONS! 1,1.) 

THEN NOPOINT:=FALSE <* CURPOS BETWEEN B POINTS *) 

ELSE 1 1 : = ! 1+1 ; <* ELSE MOVE TO NEXT SET OF POINTS *) 

END; 

(* WHEN THIS POINT IS REACHED A PROPER POINT FOR INTERPOLATION 
HAS BEEN FOUND, SO BEGIN INTERPOLATING *) 
SL:=(POSITION(. 1 1+1 ,2.)-P0SITI0N(.l 1 ,2.) 1/DPERSAM; 

XX:=CURPOS. - POSITIONC.I 1,1.); 

DXX:*XX*SL; 

1 1 1 .—1001-1; 

CORC.t 1 1 .) :*POSITION(.l 1 ,2.)+DXX; 

CURP0S:=25000. 0-1 *25.0; 

END; 

END; 

{**♦********************************************♦♦************ 

» CORSCALE....THIS IS THE FINAL SCALING PASS. THE USER IS * 

* PROMPTED FOR THE GS BEAM PATH WIDTH - FROM THE * 

* ACTUAL SITE. THE TABLE IS NORMALIZED BY THIS * 

* VALUE SO THAT IT CAN BE RESCALED BY * 

* NASA SIMULATION * 

**************************************************************) 

PROCEDURE CORSCALE; 

VAR 

1,1 1 ; INTEGER; 

OWIDTH:REAL; (* ORIGINAL LOCALIZER WIDTH *) 


BEGIN 

OWIDTH:=150.0; 

FOR I :=1 TO 1000 DO 

COR( . I . ) : =COR( . I . )/0W I DTH*0 .7; 

(* IN EXISTING NASA SOFTWARE AN ANGULAR DEVIATION OF 4-0.7 
DEGREES CAUSES A FULL SCALE DEFLECTION OF THE CO I *) 

END; 
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(1H***^***#*****«**************************************************** 

* UN I COROUT CREATES A FILE ON DISK CONTAINING THE FINISHED * 

* TABLE. 

************* 


PROCEDURE UN I COROUT; 
VAR 

1,11:1 NTEGER; 


BEGIN 

(* DEFINE FIDEFS FOR DISK 
CLOSE (OUTPUT); 

SYSTEMC'FI SYSPRINT DISK 
REWRITE(OUTPUT) ; 

FOR I :=1 TO 1000 DO 
WRITELN(C0R(.I.):10:8); 

END; 


FILE OUTPUT *) 

UC CHART C',l I); 


(*******»************ ******************* - 

* SMOOTHEND THIS PROCEDURE ’SMOOTHES' THE END OF THE TABLE * 

* SO THAT THE FIRST SMPLE WILL NOT CAUSE A BIG JUMP ON * 

* CO I NEEDLE. * 

i ****************************************************************) 


PROCEDURE SMOOTHEND; 


VAR I : INTEGER; 
SLOPE :REAL; 


BEGIN 

SLOPE:=(-COR(. 975.) )/625.0; 

FOR I :=976 TO 1000 .DO 

(* DRAW A STRAIGHT LINE FROM THE 975TH SAMPLE TO ZERO *) 
BEGIN 

COR ( . I . ) : =COR ( .975 . )+ ( I -975 ) *25*SL0PE; 

END; 


(* MAIN PROGRAM BODY 


*) 


BEGIN 

RAW IN; 

ANGLESCALE; 

DISTANCE; 

STRUCTUREOUT; 

UN I COR; 

SMOOTHEND; 

CORSCALE; 

UN I COROUT; 

END. 


it 
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PROGRAM GENPATHC I NPUT,OUTPUT) ; 

C ° NS I* THESE ARE THE CONSTANTS RELATIVE TO THE CAT I TOLLERANCE STUDY *) 
PI =3. 14593; 

NM=6076.0; 

NM2=12152.0; 

NM3= 18228.0; 

NM4=24304 . 0; 

NMH=3038.0; 

I LSPTA=24304.0; 

ILSPTB=3500.0; 

ILSPTC=1000.0; 

Z0NE2A=30.0; 

Z0NE2B a 1 5.0; 

Z0NE3B=15.0; 

ZONE3C=15.0; 


TYPE 

LABL=ARRAY(.1..13.)0F CHAR; 


(* THESE ARE THE GLOBAL VARIABLES TO BE USED *) 

VAR 

DUMMY {INTEGER; 


(* THE FOLLOWING PROCEDURE DECLARATIONS ENABLE THIS PROGRAM TO LINK TO 
THE FORTRAN PLOTTING ROUTINES IN THE LIBRARY. *) 

PROCEDURE PLOTS! IBUF:REAL;NLOC, IDEV: INTEGER); FORTRAN; 

PROCEDURE PLOT(X,Y:REAL; I PEN: INTEGER); FORTRAN; 

PROCEDURE F ACTOR (FCTR: REAL); FORTRAN; ii l xli 

PROCEDURE OFFSET(XOFF ,XFCTR,UOFF ,YFCTR:REAL ) ; FORTRAN; 

PROCEDURE WHERE(X,Y,FCTR:REAL); FORTRAN; 

PROCEDURE SYMBOL(X,Y, HEIGHT, I BCD , ANGLE s REAL; NCHAR : I NTEGER) ; FORTRAN; 

PROCEDURE NUMBER (X, Y ,HE I GHT , FNUM ,, ANGLE *EAL; NPLACE : I NTKER) ^FORTRAN ; 

PROCEDURE AXtS(X,Y:REAL; IBCO;LABL;NCHAR{ I NTEGER; SIZE, ANGLE, XMIN,OX, 

DV: REAL); FORTRAN; 

PROCEDURE LINE(XONE,YONE:REAL;NPTS, IRPT , I ALTR, I SYMs I OTEGER); FORTRAN; 

PROCEDURE SLANT(X,Y,HEIGHT:REAL; I BCD: INTEGER;ANGLE:REAL;NCHARs I NTEGER); FORTRAN, 
PROCEDURE AX; FORTRAN; 


(* THIS IS THE PROCEDURE THAT PLOTS THE AXIS AND TOLERANCE LIMITS *) 
PROCEDURE PLOTTOL; 

VAR 

1,11:1 NTEGER; 

X:REAL; 

BEGIN 

DUMMY : =0 ; 

PL0TS(DUMMY,4,1 1 ); 

AX; 

PL0TC1 .0,4.0, -3); 

PLOTCO. 0,-2. 0,2); 

PL0T(0.0,2.0,2); 

PLOTtO.O, 0.0,3); 

PL0T(8.0, 0.0,2); 

PLOT! ILSPTC/NMH,Z0NE3C/37. 5,3); 

PLOT! ILSPTB/NMH,Z0NE3B/37. 5,2); 

PLOT! ILSPTB/NMH,Z0NE2B/37. 5,3); 

PLOT! ILSPTA/NMH,Z0NE2A/37. 5,2); 

PLOT! I LSPTC/NMH ,Z0NE3C/ (-37.5) ,3); 

PLOT! ILSPTB/NMH,20NE3B/(-37.5) ,2); 

PLOT ( I LSPTB/NMH,Z0NE2B/ (-37 .5) ,3); 

PLOT! 1 LSPTA/NMH,Z0NE2A/ (-37.5) ,2) ; 

END; 

PROCEDURE CLOSEPLOT; 

BEGIN 

PLOTtO. 0,0.0,999); 

END; 
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(* THE FOLLOWING PROCEDURES GENERATE THE COURSES FOR THE THREE GENERIC 
LOCALIZER COURSES. EACH CASE PROCEDURE GENERATES A DIFFERENT COURSE 
STRUCTURE. 


PROCEDURE CASE I ; 


VAR 

X,Y,YY,A;REAL; 

(•.INTEGER; 


BEGIN 

PL0T(0. 0,0.0,35; 

FOR I :=1 TO 1000 DO 
BEGIN 


X: a ( 1-1 )*25.0; 

IF X<1 000.0 THEN Y: =-1 5. 0*C0S ( £X*P 1 5/1000.0) 

ELSE IF X<3 500.0 THEN Y:= 15*COSCCCX-1000.0)*PI )/(NM4-1000.0) ) 
ELSE Y:=( 15+(X-3500.0)*15.0/(NM4-3500.0) )* 

COS t { CX— 1 000.0)*1 .5*PI 5/CNM4-1 000.0) ); 


YY:=2.0*Y/150.0; 
WRITELNCYY: 10:8); 
PL0T(X/NMH,Y/37.5,2); 
END; 

END; 


PROCEDURE CASE 1 1; 


X,Y,YY,A:REAL; 
1,1 1 {INTEGER; 


BEGIN 

PL0T(0.0, 0.0,3); 

FOR l:»1 TO 1000 DO 
BEGIN 

X:«( 1-1 )*25.0; 

IF X < 1000.0 THEN A:=1 6.203 
ELSE IF X < 3500.0 THEN A;=15.0 

ELSE IF X < 20000.0 THEN A:*(15 + (X-3500.0)*15. 0/20804.0) 
ELSE A:=(27.229-(X-20000.0) *0.005446); 

Y s =A*COS CX/2 578.734); 

YYs=Y*2. 0/1 50.0; 

WRITELN(YY:10:8); 

P LOT (X/NMH , Y/37 . 5 ,2 ) ; 

END; 


END; 


PROCEDURE CASE I II; 
VAR 

X,Y,YY,A:REAL; 

I {INTEGER; 


BEGIN 

PLOTCO.O, 0.0,3); 

FOR I {=1 TO 1000 DO 
BEGIN 

X:=( 1-15*25.0; 

IF X<4 500.0 THEN Y:=8.493557+8.493557*C0S( (X*PI )/4500.0) 
ELSE Y;=0.0; 

YY:=Y*2. 0/1 50.0; 

WRITELN(YY : 10:8) ; 

PLOT(X/NMH, Y/37. 5,2); 

END; 

END; 

BEGIN (* MAIN PROGRAM BODY *); 

PLOTTOL; 

CASE! I ; 

CLOSEPLOT; 


END 



C THIS PROGRAM PRODUCES A PLOT FILE FOR A FILE CREATED BY 
C THE PROGRAM GSSCLV. 

CALL PLOTS (IBUF, 4,11) 

CALL AX 

CALL PLOT(1 .0,4. 0,-3) 

READ (3, 30 )Y 
X=0 

CALL PL0T(X,Y,3) 

10 F0RMAT(5X,I7) 

DO 20 1=1,999 
READ(3,30)Y 
X= I *50.0/6076.0 
Y=Y*4.0 

30 F0RMAT(F12.8) 

CALL PL0T(X,Y,2) 

20 CONTINUE 

CALL PLOTtO. 0,0.0 ,999) 

END 

SUBROUTINE AX „ c „ „ 

CALL AXI St 1. 0,2.0, 'DEGREES', 7 ,4.0,90.0, -0.5,0.25, 10.0) 

CALL AXIS(1.0,2.0, 'DISTANCE (NM) * ,-13, 7. 5, -0.0, -0.0, 0.5, 10.0) 
RETURN 
END 


COPY 41 42 43 = = = (REPLACE RECFM F LRECL 80 
FI 3 DISK 41 42 43 
FI 11 DISK PL41 42 43 
EXEC LODE GSCORPLT 



THIS PROGRAM PRODUCES A PLOT FILE FOR A FILE CREATED BY 
THE PROGRAM SCLVERT. THIS PROGRAM REQIURES THE RECORD COUNT TO 
BE CONTAINED IN THE FIRST RECORD. 

CALL PLOTS ( IBUF ,4,11) 

CALL AX 

CALL PLOT ( 1 .0 ,4 .0 ,-3 ) 

READ (3, 30 )Y 


10 


30 

20 


X=0 

CALL PL0T(X,Y,3) 
F0RMAT(5X,I7) 

DO 20 1=1 ,999 
READ(3,30)Y 
X= I *50.0/6076.0 
Y=Y*2.0 
F0RMATCF12.B) 

CALL PL0T(X,Y,2) 
CONTINUE 

CALL PLOT(0. 0,0.0, 999) 


END 

SUBROUTINE AX _ „ , „ A c , n m 

CALL AXIS(1.0,2.0, 'DEGREES' ,7,4.0,90.0,-1.0,0.5,10.0) 

CALL AXIS(1. 0,2.0, 'DISTANCE (NM) ' ,-13, 7. 5, -0.0, -0.0, 0.5, 10.0) 
RETURN 
END 


COPY 41 42 43 = = = (REPUCE RECFM F LRECL 80 
FI 3 DISK 41 42 43 
FI 1 1 DISK PL41 42 43 
EXEC LODE CORPLT 
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c * 
c * 
c * 
c * 
c * 
c * 
c * 
c * 


20 

30 

10 


C 

22 

50 

C 

23 

100 


THIS PROGRAM READS A GENERIC APPROACH PATH FRO DISK FILE 5 
AND USES THE TABLE LOOKUP METHOD TO ADD IN THE DIGITIZED 
LOCALIZER COURSE ERROR (WHICH HAS BEEN READ IN FROM DISK 
FILE 4.) THE SUBROUTINE IRSOUT IS THEN USED TO PUT THE 
OUTPUT INTO A FORM COMPATIBLE WITH THE INTELLIGENT REMOTE 
SERIAL DEVICE. VALUES OF -+2048 CORRESPOND TO FULL SCALE 
DEFLECTION OF THE LOCALIZER NEEDLE. (WHICH IS 2 DEGREES 
ACCORDING TO THE SCALING OF THE ERROR DATA ). 

DIMENSION X(1 000) ,ER( 1 000) 


DO 10 1=1 ,1000 
READ(5,20)X( I ) 
F0RMAT(3(1X,F10.4>) 

READ (4 ,30)ER( I ) 

F0RMAT(F10.8) 

CONTINUE 
XLOCER=0.0 
DO 100 1=1 ,1999 
IF (X(I).GT.25000.0) GOTO 50 
XLOCER=ER( I NT (X ( I )/25. ) ) 
XLOCER=X LOC ER* 1 024 • 0 
WRITE(2,22)XL0CER 
FORMAT ( f \ FI 0.4) 

CONTINUE 

WRITE(2,23)L0CER 
FORMAT ( * 1 , 15) 

CALL l RSOUT (XLOCER, EFLAG ) 


CONTINUE 

STOP 

END 


* 

* 

* 


* 

♦ 


* 


j 
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15 


10 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


COMMON DELT 
LOGICAL *4 BEGIN 
DELT=1 .0 
IBUF s O 

BEGIN*.TRUE. 

X*0 

Y=0 

Z=0 

CALL PL0TSCIBUF,4,11) 

CALL USRP0SC1 .0,X,Y,Z, BEGIN) 

CALL CARSC 

CALL RUNWAY 

XX=X/500.0 

YY=Y/500.0 

ZZ=Z/500.0 

CALL CARPLCXX,YY,ZZ,3) 

WRITE(6,15)X,Y,Z 

FORMATC3(1X,F10.4)) 

BEGIN*. FALSE. 

DO 10 1*1,999 
T*l 

CALL USRPOS(T,X,Y,Z, BEGIN) 

X=X/500.0 

Y=Y/500.0 

Z=Z/500.0 

CALL CARPL(X,Y,Z,2) 

WRITE(6,15)X,Y,Z 

CONTINUE 

END 


CARTESIAN PLOT 


SUBROUTINE CARPL(X,Y,Z,tPEN) 
XX=Y -0.7071 67*X+4 . 1 5 
YY=Z-0.707167*X+4 .15 
CALL PLOTCXX.YY, IPEN) 

RETURN 

END 


DRAW CARTESIAN AXIS. 


SUBROUTINE CARSC „ 

CALL AXI S(4 .15,4.15, *Y AXI S' ,6,5.85,0.0,0.0,1 .0,10.0) 
CALL AXI S(4. 15 ,4. 15, 'Z AXIS* ,6,5.85,90.0,0.0,1.0,10.0) 
CALL AXIS(4.15,4.15,*X AXIS' ,6,4.5,225.0,0.0,1.0,10.0) 


RETURN 

END 


DRAW RUNWAY 


SUBROUTINE RUNWAY 


USR00010 
USR00020 
USR00030 
USR00040 
USR00050 
USR00060 
USR00070 
USR00080 
USR00090 
USR00100 
USR001 10 
USR00120 
USR00130 
USR00140 
USR00150 
USR00160 
USR00170 
USR00180 
USR00190 
USR00200 
USR00210 
USR00220 
USR00230 
USR00240 
USR00250 
USR00260 
USR00270 
USR00280 
USR00290 
USR00300 
USR00310 
USR00320 
USR00330 
USR00340 
USR00350 
USR00360 
USR00370 
USR00380 
USR00390 
USR00400 
USR00410 
USR00420 
USR00430 
USR00440 
USR004 50 
USR00460 
USR00470 
USR00480 
USR004 90 
USR00500 
USR00510 
. USR00520 
USR00530 
USR00540 
USR00550 
USR00560 
USR00570 
USR00580 
USR00590 
USR00600 
USR00610 
USR00620 
USR00630 
USR00640 
USR00650 
USR00660 
USR00670 
USR00680 
USR00690 
USR00700 


L 


y 
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J 




10 


20 


30 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


20 . 

10 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


CALL CARPLCO. 0,0.2, 0.0,3) 
CALL CARPL(-5.0,0.2,0.0,2) 
CALL CARPL(0. 0,-0. 2,0. 0,3) 
CALL C ARP L(- 5.0, -0.2, 0.0, 2) 
CALL CARPL(0.0, -0.2, 0.0,3) 
CALL C ARP L (0.0, 0.2, 0.0, 2) 
CALL C ARP L (0.0, 0.0, 0.0, 3) 
CALL CARPL(-5.0,0.0,0.0,2) 
CALL CARPLCO .0,0. 0,0.0, 3) 
RETURN 
END 


VECTOR INTEGRATION 

THIS SUBROUTINE ACCEPTS AN ARRAY X(3,1000) AND USES 
THE TRAPAZOID METHOD TO CALCULATE THE INTEGRAL OF THE FUNCTION 
REPRESENTED IN THE ARRAY. THE INITIAL CONDITIONS FOR THE 
INTEGRATION ARE SUPPLIED THROUGH THE ARRAY I CON (3). 

DELT IS A REAL VARIABLE SPECIFYING THE SIZE OF THE INTEGRATION 
STEPS. THE CALLING FORM IS: 

SUBROUTINE VECINT(VECI N,VECOUT, I CON, DELT) 

WHERE: 

VECIN - IS THE INPUT VECTOR (3X1000) 

VECOUT - IS THE OUPUT VECTOR (3X1000) 

ICON - IS THE INITIAL CONDITIONS (3) 

CELT - IS THE DIFFERENTIAL SIZE 


REAL*4 VECIN(3,1000) 

REAL*4 VECOUT(3,1000) 

REAL *4 I CON (3) 

REAL *4 SUM 

DO 10 J=1 ,3 

SUM=ICON(J) 

DO 20 1=1 ,999 

SUM=SUM+( VEC I N( J , I )+VEC I N ( J , I +1 ) )/2*0ELT 
VECOUT(J, I )=SUM 
CONTI NUE 

VECOUT ( J , 1 000 ) =SUM+( VEC I N ( J , 999 )+VEC I N(J , 1000) )/2*DELT 
CONTINUE 
RETURN 
END 


USER POSTITION CALULATOR 


THIS ROUTINE WILL CALCULATE AN AIRCRAFTS FLIGHT PATH ACCORDING TO 
THE USERS INITIAL SPECIFICATIONS AND FLIGHT INSTRUCTIONS. 

THE ROUTINE THEN STORES (1000 SAMPLES OF) THE FLIGHT PATH IN AN 
ARRAY. THIS STEP REQUIRES AN INITIALIZING CALL TO THIS ROUTINE. 
SUBSEQUENT CALLS TO THE ROUTINE RESULTS IN THE RETURN OF THE 
POSITION AT THE TIME REQUESTED. THE CALLING FORM IS: 

SUBROUTINE USRPOS(T,X,Y,Z,BEGIN) 

WHERE* 

T - IS THE TIME (REAL) FOR WHICH THE POSITION IS TO 
BE RETURNED. 


USR00710 

USR00720 

USR00730 

USR00740 

USR00750 

USR00760 

USR00770 

USR00780 

USR00790 

USR00800 

USR0081 0 

USR00820 

USR00830 

USR00840 

USR00850 

USR00860 

USR00870 

USR00880 

USR00890 

USR00900 

USR00910 

USR00920 

USR00930 

USR00940 

USR00950 

USR00960 

USR00970 

USR00980 

USR00990 

USR01000 

USR01010 

USR01020 

USR01030 

USR01040 

USR01050 

USR01060 

USR01070 

USR01080 

USR01090 

USR01100 

USR01110 

USR01120 

USR01130 

USR01140 

USR01150 

USR01160 

USR01170 

USR01180 

USR01190 

USR01200 

USR01210 

USR01220 

USR01230 

USR01240 

USR01250 

USR01260 

US RO 1270 

USR01280 

USR01290 

USR01300 

USR01310 

USR01320 

USR01330 

USR01340 

USR01350 

USR01360 

USR01370 

USR01380 

USR01390 

USR01400 
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1 

c 


12 


C 

11 

C 

C 

20 

100 


C 

50 

110 

c 

c 

60 

C 

C 


C 


C 

1000 

c 


X - IS THE X CORRD I NATE (REAL*4 ) AT TIME T. USR01410 

Y - IS THE Y CORRD INATE (REAL *4 ) AT TIME T. USR01420 

Z - IS THE Z CORRD I NATE (REAL *4 ) AT TIME T. USR01430 

BEGIN - IS A FLAG (L0GICAL*4) INDICATING WHETHER THIS USR01440 

IS AN INITIALIZING CALL OR NOT. IF BEGIN=.TRUE.USR01450 
THE ROUTINE READS THE FLIGHT INSTRUCTIONS AND USR01460 
CALCULTES THE FIRST 1000 POINTS OF THE FLIGHT USR01470 
PATH. IF BEGIN=. FALSE. THEN THE ROUTINE SIMPLY USR01480 
ATTEMPTS TO LOOK UP IN THE PREVIOUSLY CALCULATEDUSR01490 
POSITION ARRAY THE POSITION CORRESPONDING TO USR01500 
TIME T. USR01510 

USR01520 


COMMON DELT 

DATA PI/3.141593/ ,RPD/0.017453/ ,DPR/57.29578/ 

LOG I CAL *4 BEGIN 
REAL *4 I CON (3) 

REAL*4 VEL(3 ,1000) 

REAL*4 P0S(3, 1000) 

REAL*4 TIME (1000) 

REAL *4 TIMES (1000) 

REAL *4 HEAD (1000) 

REAL *4 ACCELF(IOOO) 

REAL *4 VELF(IOOO) 

REAL *4 ACCELZ(IOOO) 

REAL *4 VELZ(IOOO) 

CONTINUE 

IF (.NOT. (BEG IN)) GOTO 1010 

DO THIS SECTION IF THIS IS AN INITIALIZATION CALL 

COUNT=0 

PS 1=0 

DO 12 I =1 ,3 
ICON( I )=0.0 
TMAX=0.0 
READ(5,11) 

SKIP HEADING CARD 
FORMAT ( 1 0X) 

READ(5,20)TlNlT,XINlT,YlNlT,ZINIT,HEADtN,VELIN 
READ ANOTHER HEADING CARD 
READ (5 ,11) 

READ INITIAL VALUES 
FORMAT(6(1X,F10.4) ) 

CONTINUE 

READ (5 ,50 ,END=1 10) (TIMES (I ) , HEAD (I ),VELF(I ),VELZ( I ) ,ACCELF ( I ) , 
♦ACCELZ(l), 1-1, 1000) 

READ AND STORE FLIGHT INSTRUCTIONS 
F0RMAT(6 (1 X,F10.4)) 

CONTINUE 

REMEMBER HOW MANY INSTRUCTIONS THERE WERE 
MCOUNT=l 

CLEAR OUT ACCELERATION ARRAY 
DO 60 J=1 ,3 
VEL(J,1000)=Q. 

CALCULATE INITIAL HEADING ANGLE 
PS I =HEAD I N*RPD*DELT 
LOAD INITIAL VELOCITIES. 

VEL( 1 , 1000 )=VELI N*COS (PS I ) 

VEL (2 , 1 000 )=VEL I N*$ I N (PS l ) 

VEL(3, 1000)=0*0 
LOAD INITIAL POSITIONS 
P0S(1 ,1000)=XI NIT 
P0S(2, 1000)=Y IN1T 
P0S(3, 1000)=ZI NIT 
INITIALIZE TIME 
TMAX=T1NIT 
CONTINUE 

CLEAR RB4AINING FLIGHT COMMAND ARRAYS 
DO 1010 I=MCOUNT, 1000 
T!MES(l)=0.0 
HEAD(l)=0.0 
VELF( I )=0.0 


USR01530 

USR01540 

USR01550 

USR01560 

US RO 1570 

USR01580 

USR01590 

USR01600 

USR01610 

USR01620 

USR01630 

USR01640 

USR01650 

USR01660 

USR01670 

USR01680 

USR01690 

USR01700 

USR01710 

USR01720 

USR01730 

USR01740 

USR01750 

USR01760 

USR01770 

US R0 1780 

USR01790 

USR01800 

USR01810 

USR01820 

USR01830 

USR01840 

USR01850 

USR01860 

USR01870 

USR01880 

USR01890 

USR01900 

USR01910 

USR01920 

USR01930 

USR01940 

USR01950 

USR01960 

USR01970 

USR01980 

USR01990 

USR02000 

USR02010 

USR02020 

USR02030 

USR02040 

USR02050 

USR02060 

USR02070 

USR02080 

USR02090 

USR02100 
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0 0-000 


010 


1100 

c 


c 

1400 


VELZC I )=0.0 
ACCELF ( 1 )=0.0 
ACCELZ( I )=0.0 

COME HERE IF THIS IS NOT AN INITIALIZATION CALL. 

CONTINUE 

IF REQUSESTED TIME IS IN THE RANGE OF ALREADY CALCULATED VALUES 
SKIP CALCUALTION PHASE 
IF (T.LT.TMAX) GOTO 2000 
CALCULATE NEW TIME INTERVAL (1000) POINTS 
DO 1500 1=1,1000 

TIME(I)=TMAX+(I-1)*Dl.~ 

TEST TO SEE IF STILL IN fL-v.'GE OF CURRENT FLIGHT 
IF (TIME(I).LT.TIMES(C0UNT+1)) GOTO 1400 
CHECH IF THERE ARE ANY MORE INSTRUCTIONS. 

IF (COUNT.GE.MCOUNT) GOTO 1400 
GET NEXT FLIGHT INSTRUCTION 
C0UNT=C0UNT+1 
CALCULATE HEADING ANGLE 
PS I =HEAD (COUNT) *DELT*RPD+PS I 

I F(FVEL.NE.VELF (COUNT) ) FVEL=FVEL+ACCELF (COUNT) 

I F(ZVEL.NE.VELZ (COUNT) ) ZVEL=ZVEL+ACCELZ (COUNT) 


USR021 10 
USR02120 
USR02130 
USR02140 
USR02150 
USR02160 
USR02170 
USR02180 
USR02190 
USR02200 
USR02210 
USR02220 
USR02230 
I NSTRCTIONUSR02240 
USR02250 
USR02260 
USR02270 
USR02280 
USR02290 
USR02300 
USR02310 
USR02320 
USR02330 


1500 

C 


2000 

C 

C 

c 


I F ( FVEL.GE.VELF (COUNT) .AND. (FVEL-ABS( ACCELF (COUNT) ) ) .LE. USR02340 
•VELF(COUNT) ) FVEL=VELF (COUNT) USR02350 

I F(ZVEL.GE.VELZ (COUNT) .AND. (ZVEL-ABS(ACCELZ(COUNT) ) ) .LE. USR02360 
*VELF (COUNT) ) ZVEL=VELF (COUNT) USR02370 

LOAD VEL ARRAY USR02380 

VEL(1,I )=VELF(COUNT)*COS(PSI) USR02390 

VEL(2,I )=VELF(COUNT)*SIN(PSI ) USR02400 

VEL(3, 1 )=VELZ(COUNT> USR02410 

CONTINUE USR02420 

INCREMENT TIME INTERVAL USR02430 

TMAX=T IME ( 1 000 )+D ELT USR02440 

LOAD INITIAL POSITION CONDITIONS USR02450 

ICONO )=P0S(1 ,1000) USR02460 

IC0N(2)=P0S(2, 1000) USR02470 

ICON(3)=POS(3,1000) USR02480 

PERFORM VECTOR INTEGRATION ON VELOCITY TO OBTAIN POS USR02490 

CALL VECI NT(VEL, POS, ICON, DELT) USR02500 

CONTINUE USR02510 

COME HERE I FT IS IN PREVIOUSLY CALCULTED RANGE USR02520 

OF POSITIONS USR02530 

LOOK UP X, Y, AND Z IN THE POSITION ARRAY USR02540 

N=l NT( (T-(TMAX-1 000.0*DELT) )/DELT) USR02550 

IF (N.EQ.O) N*1 USR02560 

X=P0S(1,N) USR02570 

Y=P0SC2,N) USR02580 

Z=P0S(3,N) USR02590 

RETURN USR02600 

DEBUG UNIT(6), TRACE USR02610 

END USR02620 
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ACONTROL ALL 

AIF .A1 EQ .? AGOTO -TELL 
ACONTROL ALL 
AX = 0 
ALOOP 5 10 
AX = AX + 1 

FI 9 DISK UCCAX CHART C 
FI 8 DISK UCGSAX CHART C 
FI 12 TAPI 

LOAD TAPETRAN (START 

AX * 0 

AY = 0 

ALOOP 8 3 

AX = AX + 1 

ALOOP 5 3 

AY - AY + 1 

FI 9 DISK GENLOCAX CHART C 
FI 8 DISK GENGSAY CHART C 
FI 12 TAPI 

LOAD TAPETRAN (START 

AY = 0 

AEXIT 

ABEGTYPE 

THIS EXEC WRITES THE TEN COURSES PREPARED FOR NASASIM TO MAG TAPE. 
THE EXEC EXPECTS TO SEE ALL TEN COURSES WITH UCC PREFIX ON THE A DISK. 
THE EXEC ALSO EXPECTS TO SEE THREE GENERIC PATH 
IT WRITES THESE TO TAPE IN ALL POSSIBLE COMBINATIONS 

AEND 
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This FORTRAN program writes the Localizer and Glide Slope error information to 
tape. 


DO 50 J=1 ,1000 
READ(9,20)LOC 
READ (8 ,20)GLS 
WRITE(12,30)L0C 
WRITE(12,30)GLS 
20 F0RMAT(F10.8) 

30 FORMAT(F10.8) 

50 CONTINUE 
STOP 
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